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Abstract 

Topological charge susceptibility Xt for pure gauge SU(3) theory at finite temperature is studied using anisotropic 
lattices. The over-improved stout-link smoothing method is utilized to calculate the topological charge. Near the phase 
transition point 'we find a rapid declining behavior for;^^; ■with values decreasing from (188(l)MeV)"^ to (67(3)MeV)"^ as 
the temperature increased from zero temperature to 1.97’^ 'which demonstrates the existence of topological excitations 
far above Tc- The 4th order cumulant C 4 of topological charge, as 'well as the ratio C 4 /xt are also investigated. Results 
of C 4 sho'w step-like behavior near Tc 'while the ratio at high temperature agrees 'with the value as predicted by the 
diluted instanton gas model. 
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1. Introduction 

Quantum Chromodynamics (QCD) is the gauge field theory that describes the rules of strong interaction among 
quarks and gluons, 'which has achieved great success in describing modern particle experiments. Higgs boson-used to 
be the last missing particle in this frame'work - 'was found in 2012 on LHC, claiming the landmark completeness of the 
Standard Model. From 2000s, lattice QCD 'which simulates QCD on finite lattices by po'werful computers, has gained 
enormous amount of results in accordance 'with experiments on accelerators, thanking for the incredible development 
of computers as 'well as the continuous progress in computational methods. Despite of the success in the spectroscopy 
of hadrons, lattice QCD researchers have made great efforts on QCD thermodynamics 'with nonzero temperature and 
chemical potential, from 'which further understanding about hadronic matters and quark-gluon plasma can be obtained 
and problems such as the physics in early universe 'with extremely high temperature and density can be addressed. 

Topology in lattice QCD has been 'widely studied for various purposes, 'while one of the most exciting applications 
is to explore the structure of QCD vacuum, 'which has not been 'well understood. The topological susceptibility Xt 
has attracted special interest since 1979, 'when the quenched 'was related to the U (1) axial anomaly and mass of 
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T]' meson through the well-known Witten-Vaneziano relation IHIIl. Furthermore, the distribution of the topological 
charge Q can be described in terms of its cumulants 0 : 


(S) 
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where e^c (0) is vacuum energy density in the 0 vacuum. The 2nd order cumulant C 2 is known as the topological 
susceptibility Xt, and C 4 - the 4th order cumulant of Q is important to lattice calculations of observables with fixed 
topo log y sector ||4l|5]. Nonzero values for C 2 „ with n > 1 indicates non-Gaussian distribution of Q ||6l|7l. From Q 


and (3 


, it is easy to see that C 4 = 




On the other hand, lattice studies of QCD at finite temperature and the phase structure of QCD have made great 
progress recently HI. In particular, the relation between chiral symmetry breaking and deconfinement transition, as 
both of them happen to be near Tc, is of great importance and has attracted many researchers. Since Xt describes 
the topological fluctuations of the vacuum in quenched situation, the behavior of Xt near Tc is expected to provide a 
further understanding on the relation between chiral symmetry breaking and confinement. 

Topological susceptibility at zero temperature for SU(3) gauge theory has been determined as (191(5)MeV)"^ 0. 
Xt at finite temperature also attracts great interests, especially around the transition temperature Tc inoHni. In 2002 , 
Gattringer et al IfT^ found the cross-over behavior of Xt in the temperature interval of O.STc ~ 1.37’^, which is in 
accordance with results of other researchers. In their work chirally improved fermion action and fermionic method 
based on Atiyah-Singer index theorem were employed to calculate the topological charge. There are basically two 
kinds of methods to extract the topological charge: fermionic and bosonic method. For bosonic method, some kind 
of smoothing procedure always be applied first to dampen the fluctuations at the UV scale, while hopefully leaving 
the long range physics unchanged. In 2010 Liischer proposed the Wilson flow method lfT4l that could be employed 
to calculate topological charge, which is renormalizable and converges to the continuum definition. More recently, 
Glaudio Bonati et al studied the dependence of 4D SU(N) gauge theories on the topological term at finite temperature 
by using cooling method iffSll . found that the 0 dependence are drastically changed across the deconfinement transition. 
Moreover, the comparison between standard cooling method and Wilson flow method were also studied in Ref. m, 
which led to equivalent results, both for average quantities and configuration by configuration. In this paper we 
present our work using over-improved stout-link smearing method uni on the anisotropic lattices with a wider range 
of temperature to provide an alternative study on this interesting subject. This smearing method is based on the 
classical instanton solution in continuum limit, which is model dependent. However, it is observed that the method 
keeps the topological charge stable and always values around the same integer obtained by the fermionic method. 
This method is also relatively cheaper than the other smoothing method. More recent studies were also carried out in 
Ref. HBHlOl and part of their results are compatible with ours. 

The diluted instanton gas model gives the ratio of C 4 and C 2 as the value -1 ED. Recent calculation at zero 
temperature gives the ratio R = -0.233(45) Q (different definition of R leads to the opposite sign with ours in that 
paper), which is in contradiction with the value predicted by the diluted instanton gas model. At high temperature 
phase(> l.lSTc) the domination of the diluted instanton gas model was confirmed through the calculation of 7>2 and 
free energy El, where 7>2 = cnl\2xf In this paper, we calculate the ratio R - CA,lxt below and above Tc with pure 
gauge configurations, trying to investigate the behavior of R around Tc using gluonic method. 

This paper is organized as follows: we first present our lattice setup and details about extracting the topological 
charge Q at finite temperature. We also discuss the over-improved stout-link method to smooth the gauge field and 
the calculation of Q by field theoretical definition employing highly improved field strength tensor ll22l . After that, 
results of Q and the susceptibility Xt^ ns well as the 4th order cumulant C 4 of topological charge and the ratio R are 
presented and discussed. Finally main results of this paper are summarized. 
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2. Lattice setting and method 


The calculation is carried out using anisotropic lattices, which have advantages of improving accuracy in lattice 
QCD for both zero and finite temperatures ll23l . The Symanzik and tadpole improvement schemes of the gauge action 
are found to have better continuum extrapolation behaviors for many physical quantities, that is, the finite lattice 
spacing effect is well suppressed by these improvements. Considering these facts, we use the following improved 
gauge action, Il24l4^ 
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where p is related to the bare QCD coupling constant go, ^ = Os/ai is the aspect ratio for anisotropy (we take ^ - 5 
in this work), and m, are the tadpole improvement parameters of spatial and temporal gauge links respectively. 
Qc = YjC |ReTr(l - Wc), with 14c referring to the path ordered product of link variables along a closed contour C 
on the lattice. Qsp includes the sum over all spatial plaquettes on the lattice, Dtp includes the temporal plaquettes, Qsr 
includes the product of link variables about planar 2 x 1 spatial rectangular loops, and Qstr refers to the short temporal 
rectangles (one temporal link, two spatial). Practically, u, is set to 1, and Us is defined by the expectation value of the 
spatial plaquette P^s'- Us - ^ . Besides, it has been shown in Ref. IIZTl that the renormalization effects of the 

anisotropy ^ is ignorable as p varies. We also have used Wilson flow method to calculate the renormalized anisotropy 
according to Ref. Il28l at different temperature with p = 3.2, and the results show that the differences between the 
renormalized ^ and bare ones are always less than 2%. We adopt anisotropic lattices with smaller a, compared to as, 
so that we can investigate higher temperature with small spatial lattice size. 


2.1. The definition of temperature 

The temperature on lattice is defined as follow: 


1 

T = -, 
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where N, is the temporal lattice size. T can be changed by varying either Nt or the coupling constant p, which is 
related directly to the lattice spacing a,. The critical temperature is determined for a given lattice Nt - 24 after the 
critical coupling Pc has been determined. The order parameter for determining Pc is chosen to be the susceptibility 
of the Polyakov line, which is defined as 


T. - <©') 

where 0 is the Z 3 rotated Polyakov line defined via: 
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and P represents the trace of the spatially averaged Polyakov line for each gauge configuration. After a rough scan, a 
more refined study for the peak position of the susceptibility gives the critical coupling constant Pc - 2.808, which 
corresponds to the critical temperature Tc ~ 0.724rg' = 296MeV Ell. Here ro is the hadronic scale parameter and 
we take = 410(20)MeV. 

Considering both finite volume effects and good resolution of temporal lattice at T ~ 2Tc, we set p - 3.2 in the 
study of topological susceptibility at finite temperature. The corresponding lattice spacing as is obtained by calculating 
the static quark potential V (r) on an anisotropic lattice 24^ x 128. The fitting result of string tension cr reads: 
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taking r^^ =410 (20)MeV, we have as — 0.0878 (4) fm. Comparing to as at Pc - 2.808 and A, = 24, the approximate 
values for A, that corresponding to and 2Tc at fixed p = 3.2 can be calculated and the results are A, ~ 38 and 
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~ 19 respectively. More parameter details can be found in Ref. UTij . It should be pointed out that the difference 
between TcijS - 2.808) and TciP - 3.2) is negligible due to the application of the improved gauge action. 

In this study we change the temperature by varying temporal lattice size Nf Working on the same /? which 
corresponding to the same a, makes all our lattices share the same spatial volume while the temperature T varies in 
a wide range. This avoids the influence from possible finite volume effects arising from different spatial volumes. 
Setting p - 3.2, we generated a series of lattice 24^ x Nt with Nt =20, 24, 28, 32, 36, 40, 44, 48, 60, 80 and 128, 
which cover the range of T e [0.37’c, l.9Tc\, and the spatial size Lj ^ 2.1fm. For each anisotropic lattice with a 
fixed value of Nt, we sampled roughly 1000 configurations(1900 for Nt - 60), after 10000 sweeps from cold start and 
with 500 sweeps between each sampling. Here one sweep consists of a composite update of 1 pseudo heat-bath and 
5 over-relaxation procedures over all link variables. 5000 samples of 24^ isotropic lattice with the same were also 
generated and measured as the results at zero temperature. We roughly estimated the integrated autocorrelation time 
of topological charge tq, which is about 200 sweeps at zero temperature and about 440 sweeps at \.9Tc. 


2.2. Over-improved stout-link method 

It is well-known that the topological charge Q calculated directly from a typical lattice conhguration by the held 
theoretic dehnition ll29l : 

J (9) 

j ? ^fiv — 'i^^iivpa'^pa-, (10) 

will not be an integer in general, which is obviously conhicting with the continuum situation. Thus, some cooling 
or smearing procedure is essentially to be taken on original conhgurations before calculating Q by bosonic methods, 
aiming to suppress the ultraviolet Huctuations. Highly improved lattice held-strength tensor ll22l is also helpful to 
obtain an integer Q on the lattice, as we used 1 x 1, 2 x 2 and 3x3 3-loop improved which is practically good 
enough. 

All cooling or smearing methods are based on an approximation to the continuum gauge held action: 



and the lattice version of Sg contains combination of plaquette and larger Wilson loops. However, the discretization 
error always exists, which is harmful to the topological objects and smoothing procedure. So in these methods the 
destruction of topological structure is unavoidable. For the purpose of the study on topology susceptibility, we adopt 
the over-improved stout-link smearing method developed by Moran and Leinweber Oil to hlter the lattice, which is 
proved to be efficient in preserving instanton-like objects. 

Over-improved stout-link smearing method introduces the over-improved parameter e ll3Tll into stout-link smear¬ 
ing algorithm Il24l . The improved action with instanton solution substituted reads: 


5'''"'(e) = 
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where pi„st is the instanton size. Negative e leads to a positive leading order error term of p'which means instanton- 
like objects will be preserved when the action is reduced in smoothing procedure. The modified link combination used 
by over-improved stout-link method is as follow: 

Cp (x) = Psm ^ I (1 1 touching Up (x)) 

—j^(l X 2 -H 2 X 1 paths touching Up (.r))|. (13) 

There are three free parameters in the over-improved stout-link smearing method: the over-improved parameter 
e, the smearing parameter psm, and Wsm - the number of smoothing steps taken on each configuration. Moran and 
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Figure 1: Smooth tests for different p aiN, - 40. Tests at other temperatures show similar situation, except that at 
lower temperatures, more steps are needed to extract the nearly integer values of Q. 


Leinweber had practically determined e = -0.25 to maximize the life of instantons under iterative smearing procedure, 
which we also kept hxed in our study. We check this on several configurations at different temperature and conhrmed 
that topological charge approaches integral number quickly as long as Wsm ^ 10 ~ 40 in most situation, and keeps 
almost the same value for thousands of steps, which will be presented at the last of this sector. 

For the other two parameters psm and we perform some tests to find suitable values for them. Five config¬ 
urations are investigated for each temperature and smoothed 200 times for psm = 0.01 ~ 0.08, and the results of 
topological charge for some some typical values of psm are shown in Figure It’s obvious that the larger psm be¬ 
comes, the faster the topological charge becomes an integer. However, in practice we found that when ps^, is equal to 
or larger than about 0.08, smoothing leads to completely unstable results, as shown in lower right panel in Figure [1] 
which is the same situation for all Nt. Finally we set Wsm = 40 and psm = 0.05 for all lattices (except for Nt - 128 we 
set Psm = 0.07, keeping Wsm = 40 ). 




Figure 2: Stability tests at T = l.l9Tc (left) and T - 0.19Tc (right) on 10 conhgurations with nonzero Q. 

By the way, we check the stability of the over-improved stout-link smearing method, which is supposed to be an 
effective and cheaper method (compared to the fermionic method) to suppressing the ultraviolet fluctuations on the 
lattices while preserving topological structures as long-lived as possible during smoothing. One of the most important 
problems of smoothing-type methods is the instability of topological charge when smoothing is carried out. We made 
several tests to check this, such as 10000 steps smearing on different configurations with nonzero topological charge 
Q both below and above T^. Part of the results are shown in Figure It indicates that, within the over-improved 
stout-link smearing method, with proper parameters and improved field strength tensor adopted, one can preserved 
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instanton-like objects for long enough time without causing the system to fall into the trivial topological sector. 
Therefore, the result for the topological charge is robust and almost independent of the number of smoothing steps 
once an integral value for Q is obtained for most of conhgurations. Notice that we measure the topological charge 
after every 500 steps of smoothing in stability tests while for most of conhgurations Q is almost an integer after about 
40 steps. 

3. Results and discussion 

We extract 1000 topological charge data for each temperature on anisotropic lattices and 5000 samples on 24“* 
isotropic lattice, which are all obtained after 40 steps of smoothing procedures. Most of these values are nearly 
integers. Those Q values that deviate larger than 0.1 from an integer - which we called "bad points" - are less 
than 4% for each lattice setting. The number of "bad points" grows slightly as the temperature approaches zero and 
for Nt - 128 (corresponding to nearly zero temperature) we set psm = 0.07 to ensure that 40 steps of smooth are 
practically enough to avoid too much "bad points". 





Figure 3; Distribution of the topological charge at different temperature. 

Firstly, we display the histograms of the topological charge Q for different temperatures in Figure As the 
temperature increases, distribution of Q becomes narrower as expected, that the topological excitations are suppressed 
above Tc- We ht the histograms of Q with a Gaussian function which are denoted as the blue lines in Figure]^ The 
deviation from the Gaussian distribution comes from nonzero higher order cumulants and limited sampling. 

Next, the susceptibility Xt the 4th-order cumulant C 4 of Q as well as the ratio R are worked out at different 
temperatures, and the results are shown in Figure]^ Figure|^and Table[2where the errors are statistical and estimated 
by jackknife method. These quantities from original topological charge data as well as Q after rounding off are 
calculated, which make almost no differences. So the results presented are based on original Q measured. 

The susceptibility at T = 0 reads = (188(l)MeV)^ which is consistent with former results by other methods. It 
also clearly shows that topological excitation exists even aiT - l.9Tc, where= (67MeV)"^ and 4.2% conhgurations 
have nonzero topological charge. That means instanton-like structures may exist even above Tc to nearly 2Tc support¬ 
ing the conclusion of Gattringer ef a/. As Figure|^shows, a crossover behavior fordecreasing from (188(1 )MeV)"* to 
(67(3)MeV)'* is found. For all conhgurations below Tc (from 0 to 0.95Tc),Xt stays around (177MeV)"^ ~ (189MeV)"^. 
From 24^ x 36 conhgurations corresponding to l.QSTc, Xt starts to decline as the temperature increases. Notice that 
Figure|^shows the value of;t''^"^, and;^^, decreases much faster than seen above Tc, consistent with results in Ref. ina. 
For example, at T = 1.057’^ where = (148MeV)"^ which is 38.4% of its zero-temperature value, while at T = 1.97’c, 
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Xt - (67MeV)'^ has dropped to about only 1.6% of its zero-temperature value. The value of;^^, remains almost constant 
below Tc and declines rather rapidly from 0.95rc to l.QSTc signaling for the phase transition. 



Figure 4: Results for;^f’^"^ vs. temperature 



Figure 5; of Q. For convenient, we show |c 4 |*^"^ with the same sign of 


C4 


As shown in Table[2 R approaches -1 from below in the interval [1.057’c,1.197'c] while above 1.35Tc R values not 
far away from -1, in agreement with the result of Ref. 021 . in which b 2 - Rf 12 approaches -1/12 from below as T 
goes beyond or so. This is also consistent with the value predicted by the diluted instanton gas model. Similar 

coincidence has also been found in a recent separate lattice study 0?l . Besides, although the error bars grow larger 
as the temperature approaches zero, C 4 shows a step-like behavior near Tc that could be seen in Figure However, 
these behaviors couldn’t be the proof of the existence of discontinuity between O.ShTc and 0.957’^, that may be related 
to chiral and deconfinement phase transition, and should be carefully examined with more accurate calculations and 
finer temperature resolution. It is also found that, for pure gauge theories, the calculations of C 4 and R do need much 
more samples to suppress the static errors, so we can’t make any conclusion from the results of C 4 and R near 7 = 0. 
The statistical errors of Xt C 4 could be evaluated according to Q, which indeed shows much noisier signals of C 4 
than those of Xt near zero temperature. 
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config. 

r/MeV 

TIT, 

xr/Mev 

C4/MeV4 

R(c4lx) 

24^ X 20 

562 

1.90 

67(3) 

-(69(6))" 

-1.09(26) 

24^ X 24 

468 

1.58 

82(2) 

-(76(3))" 

-0.75(10) 

24^ X 28 

401 

1.35 

98(2) 

-(93(4))" 

-0.81(11) 

24^ X 32 

351 

1.19 

125(2) 

-(141(10))" 

-1.58(36) 

24^ X 36 

312 

1.05 

148(2) 

-(180(14))" 

-2.20(60) 

24^ X 40 

281 

0.95 

182(2) 

-(198(33))" 

-1.38(84) 

24^ X 44 

255 

0.86 

187(2) 

(189(29))" 

1.05(62) 

24^ X 48 

234 

0.79 

187(2) 

(184(33))" 

0.94(68) 

24^ X 60 

187 

0.63 

184(1) 

(175(37))" 

0.83(72) 

24^ X 80 

140 

0.47 

183(2) 

(138(121))" 

0.31(1.16) 

24^ X 128 

88 

0.30 

189(2) 

(145(246))" 

0.32(2.78) 

24^ X 24 

0 

0 

188(1) 

-(223(33))" 

-2.00(1.14) 


Table 1; ^ ^4 and R at different temperatures. Last row is corresponding to isotropic configurations at T = 0. 


4. Summary 

Topological susceptibility and the 4th order cumulant of topological charge in pure SU(3) gauge theory are cal¬ 
culated using anisotropic lattices from zero temperature to 1.97’c. The computation utilizes over-improved stout-link 
smearing method, employing 3-loop highly improved field strength tensor. The results obtained are consistent with 
than existing results in literatures using other methods, while extended to higher temperature. Once appropriate pa¬ 
rameters are adopted, it is confirmed that the value of the topological charge is stable and robust during the smoothing 
procedure after becoming an integer. We also confirm the existence of topological excitations around \.9Tc. The 
asymptotic behavior of the ratio c/[lxt above 1.19rc is consistent with the diluted instanton gas model. Furthermore, 
a step-like behavior of C 4 around inspires further investigation of high order cumulants of topological charge. 
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